% 创建时间：2025-06-22 12:26:39
% 作者：兪馳
clc;
clear;
A=[-7,2,0,-1;
    1,-4,2,1;
    2,-1,-6,-1;
    -1,-1,0,4];

syms k,
A=sym(A);
F=mpowersym(A,k)
F=collect(F,(-6)^k)

% disp('-------------------------------------------')
% simplify(A^123456-subs(F,k,12345))
% syms x;
% A3=funmsym(sym(A),sqrt(x),x),
% simplify(A3-subs(F,k,1/2))
% disp('-------------------------------------------')
% syms k x;
% F1=funmsym(A,x^k,x),
% simplify(F-F1)
